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Abstract: A new technique for improving the precision of measurements of SUSY par- 
ticle masses at the LHC is introduced. The technique involves kinematic fitting of events 
with two fully identified decay chains. We incorporate both event E™ tss constraints and 
independent constraints provided by kinematic end-points in experiment invariant mass 
distributions of SUSY decay products. Incorporation of the event specific information 
maximises the information used in the fit and is shown to reduce the mass measurement 
uncertainites by ~ 30% compared to conventional fitting of experiment end-point con- 
straints for the SPSla benchmark model. 
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1. Introduction 

Measurement of SUSY particle masses in R-parity conserving SUSY events at hadron 
colliders such as the LHC is complicated by a number of factors. First, imperfect detector 
hermeticity close to the beam-pipe and our ignorance of the initial state parton momentum 
prevent measurement of the final state momentum in the z(beam)-direction. This removes 
one kinematic constraint on the masses of particles generated in any given event. Second, 
the presence of an invisible Lightest SUSY Particle (LSP) at the end of each decay chain 
generates eight unknown four- momentum components in every event. As a result the 
kinematics of SUSY events at hadron colliders are typically highly under-constrained. 

Several approaches to this general problem have been documented. Given a sufficiently 
long decay chain, constraints on analytical combinations of SUSY particle ('sparticle') 
masses can be obtained from the positions of end-points in distributions of invariant masses 
of combinations of visible SUSY-decay products (jets, leptons etc.) IJ. The system of 
equations may be solved with a numerical fit to obtain the individual masses, if enough 
constraints are provided ||, g [3| . 

An alternative class of techniques which does not rely on the end-point information is 
also available. These techniques vary the assumed LSP momenta in a given event subject 
to event kinematic constraints in order to construct variables sensitive to the masses of 
sparticles present in the event. One such technique involves the construction of the Mt2 
variable [Q, [2], 0], which is usually defined for events where the same decay chain appears in 
both 'legs' of the event. For given assumed LSP transverse momenta the maximum value 
of the transverse masses of the two legs is calculated. This quantity is then minimised 
by varying the LSP momenta subject to the event E™ lss constraints to give the value of 
Mt2- The end-point of the Mt2 distribution depends on the test LSP mass M tes t and has 
recently been shown to display a kink structure when Mt es t is equal to the true LSP mass 
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B H This shows that E™ lss constraints are important for obtaining the absolute 
mass scale of SUSY particles. 

Another, more involved technique of this kind is the mass relation method [10]. In 
this method the masses of sparticles in the decay chain are calculated for a range of LSP 
momenta. When the decay chain present in the event involves a large number of SUSY 
states, the calculated sparticle masses are constrained to lie on a hypersurface in the mass 
parameter space, and the intersection of several such hypersurfaces from different events 
determines the true sparticle masses. Alternatively the existence of solutions to the system 
of mass-shell conditions provided by under-constrained events can be used to constrain the 
individual masses [pd|| . 

The techniques described above all seek to determine the masses of particles appearing 
in kinematically under-constrained events by combining information from multiple events, 
either explicitly, as in the mass-relation method, or implicitly by obtaining constraints from 
distributions of quantities derived from individual events. This paper describes a simple 
extension of these techniques (the hybrid method) in which information from multiple events 
is combined with event information in a kinematic fit to fully reconstruct events and hence 
constrain the masses of individual SUSY particles. In this paper, we especially emphasize 
a case where both legs contain the same cascade decay chain, so that events are more 
prominent over the background, and fewer unknown mass parameters are required. It is also 
possible to include E™ tss constraints for this case, which provide independent information 
on the LSP mass. Full reconstruction of SUSY particle kinematics with this technique could 
also be useful for other purposes such as measuring the spin-statistics of SUSY states. 

Section § introduces the new technique while Section ||| discusses a particular example 
of its application and illustrates the improvement in mass measurement precision obtained 
for a specific benchmark SUSY model. Section ^| concludes and outlines directions for 
future work. 



2. Description of technique 

To simplify the description, we first examine the case where a specific decay chain appears 
in both legs of the event. Such 'symmetric' events will possess a considerably smaller 
branching ratio than 'single-leg' events containing only one such decay chain but the addi- 
tional decay products can lead to striking signatures which strongly suppress the SM and 
SUSY backgrounds. Conversely, if such events are observable then it is also probable that 
significant numbers of single-leg events can also be observed. Both classes of event are used 
by the technique described here. 

Our strategy is as follows. We start with constraints on sparticle masses obtained from 
kinematic end-points in distributions of invariant mass combinations obtained from single- 
leg events, as described in Section |l|. Because one of each of these distributions is obtained 
from the experimental data-set considered we call these constraints 'experiment- wise'. In 
addition to such constraints, we possess 'event-wise' constraints which are obtained from 
mass-shell conditions provided by the event visible decay products and LSP momenta, 
where the latter are further constrained by the measured event E™ lss components. A 
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kinematic fit to both the experiment- wise and event-wise constraints can possess fewer 
degrees-of-freedom than the simple experiment-wise end-point fit, potentially improving 
the sparticle mass precision. 

In practice event visible momenta and EJp tss values may be measured with limited pre- 
cision causing each individual kinematic fit to generate mass values which potentially devi- 
ate more from the true masses than the equivalent quantities obtained from the experiment- 
wise end-point fit. Nevertheless because the data-set provides us with a number of uncor- 
rected candidate events, the means of all the fitted event mass values can indeed improve 
the mass measurement precision. The situation is analagous to measuring the mass of a 
resonance decaying to visible decay products - the mass value obtained from one event 
may be relatively inaccurate but the uncertainty is reduced by considering the mean mass 
obtained from all events. 

It is instructive to consider at this point the roles played by the different constraints 
in the event- wise fits in this simplified case of symmetric events. Each leg of the event 
contributes four unknown LSP four-momentum components to the degrees-of-freedom of 
the event, while for a decay chain consisting of n steps each leg contributes n mass-shell 
conditions. Let us define the number of degrees-of-freedom of the experiment- wise end- 
point fit to be d en d = n mass — c en d, where n mass (= n) is the number of sparticle masses and 
c end is the number of uncorrelated end-point constraints. In this case the extra number of 
degrees-of-freedom generated by including the event -wise information is d ev t — fi mom — c ev t , 
where n mom is the number of unknown LSP four-momentum components and c evt is the 
number of event-wise constraints. The condition for inclusion of the event-wise information 
to potentially improve the mass precision is d ev t < 0. 

For symmetric events d ev t is given by 4 — n for one leg and respectively 6 — 2n or 
8 — 2n for two legs with or without E™ lss constraints applied. Consequently if n > 4 even 
considering one leg alone in conjunction with the end-point constraints potentially improves 
the mass measurement precisions, while if n > 3 improved precisions can potentially be 
obtained by using both legs with E™ lss constraints. 

It should be noted that given an unambiguous assignment of visible final states to decay 
chain steps, the technique described here could be extended to non-symmetric events where 
some/all of the SUSY states appear in both legs of the event. In this case, for rt\ and ri2 
steps in legs 1 and 2, d ev t is given by 6 — (n\ + 712) or 8 — (n\ + 712) when including or 
excluding the event- wise E™ lss constraints. Consequently if n\ + ri2 > 6 improved mass 
measurement precisions can be obtained by using event-wise fits with E™ lss constraints 
applied. 

3. Example: n = 4 step symmetric decay chain 
3.1 General discussion 

Let us now consider an example of mass measurement using symmetric events, for the 
specific case of n = 4. The decay chain present in both legs is 

5 — > 7c — ► (5bc — ► aabc, (3-1) 
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where Greek letters denote SUSY states, Roman letters denote visible SM decay products 
and a is the LSP. Denoting the two legs of the event with subscripts 1 and 2, the eight 
mass-shell conditions are: 

(p(ai) +p(b 1 ) +p(ci) +p{ai)f = (p(a 2 ) + p{b 2 ) + p(c 2 ) + p(a 2 )) 2 = m$, 

(p(ai) +p(bi) +p{ai)) 2 = (p(o 2 ) +p{b2) + p(a 2 )) 2 = m*, 

(p(oi) + p(a 1 )) 2 = (p{a 2 ) + p(a 2 )) 2 = nip, 

{p( ai )) 2 = (p(a 2 )) 2 = m 2 a . (3.2) 

To these constraints should be added the two constraints provided by the event E™ tss 
components: 

p x (ai)+Px(a 2 ) = E™ ss , 

p y (a 1 )+ Py (a 2 )=E™ iss , (3.3) 

thus giving ten event- wise constraints in total (i.e. c ev t = 10). The number of unknown 
parameters appearing in these constraints is twelve - eight LSP four-momentum compo- 
nents together with four unknown sparticle masses. With the definitions from Section || the 
number of extra unknown parameters present in the event-wise fit relative to the end-point 
fit is therefore n mom = 8. Consequently d evt = —2 indicating a potential gain over the 
end-point fit. In practice the four mass-shell conditions for each leg may be solved analyti- 
cally to give a locally invertible map from the four sparticle masses to the four components 
of the LSP four momentum. In this case c ev t = 2 and n mom = 0, however d ev t is clearly 
unchanged. 

It is also interesting to consider the above example when EJp iss constraints are not 
included in the event-wise fit. In this case n mom = 8 and c ev t = 8 (or equivalently n mom = 
and c ev t = if solving the mass-shell conditions) and hence d ev t = 0. Each leg of the 
event is independent and the event-wise constraints just map sparticle masses to LSP four 
momenta. The event-wise fit is therefore equivalent to solving the kinematic end-point 
mass constraints and no gain in mass precision is obtained. The situation changes when 
the EJp iss constraints are applied because in this case the kinematics of the two legs are 
connected. 

Note that an alternative, related, approach to this n = 4 problem can also be taken. 
Namely, one merely utilises the equivalence of the fitted masses in the two legs of the event. 
In this case the constraints in Eqns. |3.2| and |3.3| reduce in number to six (since the masses 
appearing in Eqn. [T^ are a priori unconstrained) and hence two degrees-of-freedom remain 
(d ev t = 2). Thus each event is under-constrained. The LSP four-momenta which satisfy 
the constraints for a given event define a two-dimensional surface in the four-dimensional 
m a, m /3, m 'y, space. If a few events of this kind can be found, it is possible to solve 
for the masses, which are given by the coordinates of the point of intersection of all the 
event 2D surfaces in the 4D space. This approach is similar to the mass relation method 
described in the introduction and will be studied in full in a future paper. 
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3.2 Concrete example: decays at the SPSla mSUGRA benchmark point 

Let us now examine in detail a concrete realistion of the decay chain discussed above. At 
mSUGRA point SPSla there is a significant branching ratio for the decay chain 

q L - X2Q - hk -> (3.4) 

This chain provides 5 kinematic end-point mass constraints from invariant mass combina- 
tions of jets and leptons [12|: 

• m(ll) max = 77.08 ± 0.08(scale) ± 0.05(stat) GeV 

• m(llq) max = 431.1 ± 4.3(scale) ± 2.4(stat) GeV 

• m(llq) min = 203.0 ± 2.0(scale) ± 2.8(stat) GeV 

• m{lq)™ x = 380.3 ± 3.8(scale) ± 1.8(stat) GeV 

• m(lq)™ ax = 302.1 ± 3.0(scale) ± 1.5(stat) GeV 

For this study unbiased samples equivalent to 100 fb _1 (one Monte Carlo 'experiment') 



of SPSla signal events and ti background events were generated with HERWIG 6.4 14] 
and passed to a generic LHC detector simulation [15]. A lepton reconstruction efficiency 
of 90% was assumed. 

Events were selected in which the above decay chain appears in both legs of the event 
with the following requirements: 

• N jet > 2, with prO'2) > 100 GeV, 

. M e//2 = E™ iss +pr(il) +Pt0'2) > 100 GeV, 

• E^ iss > max(100 GeV,0.2M e//2 ), 

• Ni ep = 4, where lep = e/^(isolated) and pt(W) > 6 GeV, 

• 2 Opposite Sign Same Flavour (OSSF) lepton pairs. If the pairs are of different 
flavour both pairs must have mill) < m(ll) max . If both pairs are of the same flavour 
then one and only one of the two possible pairings must give two mill) values which 
are both less than m(ll) max . These pairings allocate the leptons to each leg of the 
event. 

• One and only one possible pairing of the two leading jets with the two OSSF lepton 
pairs must give two m(llq) values less than m(llq) max . These pairings allocate the 
jets to each leg of the event. 

• For each inferred leg of the event the maximum(minimum) of the two m{lq) values 
must be less than m (^)^o)- This ordering allocates the leptons to the near and 
far P] positions in the decay chain. 
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The requirement of 4-leptons in two OSSF pairs and two high-p-r jets consistent with 
kinematic end-points, together with large E™ lss , is effective at removing the majority of 
SM and SUSY backgrounds (see below). 

Each selected event was fitted with MINUIT ]16| . Free parameters were taken to be 
the four masses appearing in the decay chain: m(qx), rn{x^)i m(lji) and m(xi)- In the 



spirit of the discussion in Section 3.1 the mass-shell conditions and measured momenta of 
the visible decay products for each leg were solved to determine the LSP four-momenta, 
giving two solutions for each leg. The x 2 minimisation function was defined by: 



x 2 



'm(ll)? v r-m(ll)™r 2 



&m(ll) max 

rn(Uq)™*-m(llq)%%r) 2 | ( m(llq)™ n - mjllq)™ } 2 



a m(llq) max J \ a m(llq) min 



+ f rn(lq)^-m(lq)^ \ 2 | U(lq)^ - mgggg pt \ 2 

+ ( PAxjjl)) + P*(xg(2)) - Er ss Y | / Py(xg(l)) + Py(xg(2)) - E™» \ 2 ^ (3 5) 

y (T^miss / y a E rniss J 

where evt denotes an expected end-point value derived from the masses in the event-wise 
fit with the formulae of Ref. ||, and expt denotes a 'measured' experiment-wise end-point 
value. The uncertainties a in these 'measured' endpoints were those quoted above. The 
uncertainties on the measurements of the x and y components of E™ lss , a E m,iss and a E m.iss, 
were given by 0.5^/£^ um where E^ m is the scalar sum of jet pr of the event. This function 
incorporating both event-wise EJp tss constraints and experiment-wise end-point constraints 
was evaluated for each of the four pairs of Xx momentum solutions obtained from solving 
the leg mass-shell conditions. Fitted masses were obtained when x 2 was minimsed for the 
event. Fitted masses were used in the subsequent analysis only if MINUIT judged the fit 
to have converged and Xmin < 35.0. 

Following application of the selection cuts described above and the requirements of fit 
convergence and low fit Xmin ^ SUSY 'signal' events with the above decay chain appearing 
in both legs were observed. 4 SUSY background events were observed, consisting of the 
above decay chain in both legs but with one or two leptonically decaying staus produced in 
the decays of the X2 S - N° tt background events were observed in 100 fb _1 equivalent data. 
More SM background events may be expected in a real experiment, given that effects such 
as charge and lepton mis-identification are not included in the fast detector simulation. 
The use of full GEANT detector simulation is required to model correctly these effects, which 
is beyond the scope of this paper, however given the hard kinematic selection cuts it is 
reasonable to assume that they are at least smaller than the negligible ti background 
considered above. 

Each event-wise fit generated one set of values for the sparticle masses, namely those 
values which minimise the broad x 2 function in Eqn. |3.5| . The distributions of these values 
for one Monte Carlo experiment are shown in Fig. |l[ The widths of the distributions are 
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Figure 1: Distributions of sparticle masses obtained from event- wise fits, for one MC experiment. 
Each entry is obtained by minimising the \ 2 function shown in Eqn. 3.5 



greater than those which would be obtained were the EJp %ss constraints excluded from the 
fits - in that case each event-wise fit becomes equivalent to the experiment-wise end-point 
fit as explained in Section [Tl] and each distribution becomes a delta-function located at the 
equivalent mass value. The key point is that despite the fact that the mass distributions in 
Fig. H are broader when E™ tss constraints are used, the means of the distributions measure 
the sparticle masses more accurately. This may seem counter-intuitive, however the reason 
is clear. The EJp iss components for each event are measured with only limited precision, 
causing a potential shift in the minimum of the event x 2 function away from the true 
sparticle masses. There is no way to correct for this effect on an event-wise basis. By 
correctly parameterising the E™ ss measurement precision in Eqn. 3J3 however we ensure 
that such shifts are on average unbiased when considering all events in the experiment and 
hence the means of the mass distributions compensate for this resolution effect, while also 
making use of the reduction in degrees-of-freedom provided by the EJp tss constraints. This 
compensation for resolution effects at the level of the experiment rather than individual 
events is also implicit in the end-point method where gaussian-smeared end-point functions 
are used when fitting the end-points in order to compensate for experimental jet and lepton 
energy resolutions. 

In order to demonstrate the performance of the technique and judge the uncertainties 
in the measurements the above procedure was repeated for 100 Monte Carlo experiments. 
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Figure 2: Likelihood distributions of sparticle masses obtained from 100 MC experiments. Each 
entry is the mean of an experiment-wise mass histogram such as those in Fig. [l]. 



For each experiment, kinematic end-point positions were sampled from gaussians with 
means and sigmas given by the means and uncertainties listed above. The five sampled 
end-point positions for each experiment were solved simultaneously with a MINUIT fit 
to give initial mass values for input to the MINUIT event-wise kinematic fits. For each 
experiment relative jet(lepton) energy scale values were sampled from gaussians of width 
1%(0.1%) reflecting likely ultimate energy scale uncertainties at the LHC. Each experiment 
generated a set of sparticle mass histograms similar to those shown in Fig. [I]. The means 
of these histograms for the 100 MC experiments were then used to construct likelihood 
histograms for the masses, shown in Fig. The standard deviations of these histograms 
were taken to provide the uncertainties on the sparticle mass measurements. 

Unbiased MC data equivalent to only one 100 fb _1 experiment were available for this 
study. For this reason the same events were used for each MC experiment, with just 
the end-point values and jet/lepton energy scales varying. The additional uncertainties 
in the final mass values expected from varying event samples were estimated from the 
mean statistical uncertainties in the mean experiment mass values as extracted from the 
event-wise distributions such as those shown in Fig. |]. We evaluated the experiment-by- 
experiment spread due to varying event samples as a/y/n, where a is the RMS of the 
event-wise distributions as shown in Fig. [I], and n is the number of entries in each plot. 
These additional contributions were added in quadrature to the uncertainties obtained 
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State 


Input 


End-Point Fit 


Hybrid Method, E™ iss 


Hybrid Method, no E™ iss 


Mean 


Error 


Mean 


Error 


Mean 


Error 


X? 


96.05 


96.5 


8.0 


95.8(92.2) 


5.3(5.5) 


97.7(96.9) 


7.6(8.0) 


Ir 


142.97 


143.3 


7.9 


142.2(138.7) 


5.4(5.6) 


144.5(143.8) 


7.8(8.1) 


xl 


176.81 


177.2 


7.7 


176.4(172.8) 


5.3(5.4) 


178.4(177.6) 


7.6(7.9) 




537.2-543.0 


540.4 


12.6 


540.7(534.8) 


8.5(8.7) 


542.9(541.4) 


12.2(12.7) 



Table 1: Summary of mass measurement precisions for SPSla states. Column 2 lists masses used 
in the HERWIG generator, Columns 3 and 4 the fitted masses and uncertainties obtained from the 
conventional fit to kinematic end-points, Columns 5 and 6 the equivalent values obtained with 
the new technique and Columns 7 and 8 the equivalent values obtained with the new technique 
excluding E™ zss constraints. Figures in parentheses are those obtained with the biased sample of 
non-repeated events. All masses are in GeV. The quoted mass range for qx excludes b squarks, 
which are produced less readily than the light squarks. 



from the study. This approximation was checked with a second sample of SPSla events 
equivalent to 100 different MC experiments, biased to force gluinos to decay to qi, b or t, 
q~L to decay to x® an d X2 to decay to e or /}. Excluding the SUSY background from this 
sample was estimated to bias the mean mass values by < 1% and increase the uncertainties 
by ^5%. The uncertainties increase when the SUSY background is excluded because the 
effect of the decrease in event statistics outweighs the reduction in bias caused by excluding 
the primarily f background events, which have similar kinematics to the e and jl signal 
events. 

The results of this study are summarised in Table |lj. For comparison purposes the 
analysis was initially carried out with the EJp tss constraints removed from the x 2 function. 
The measurement precisions are consistent with those obtained from the conventional end- 
point fitting method, as expected following the reasoning outlined in Section 3.1. The 
analysis was then repeated including the E™ lss constraints, giving an overall improvement 
in sparticle mass precisions ~ 30% for all four masses considered. A similar improvement 
was found when using the biased sample of non-repeated events for different experiments. 
The measurement of mass differences is also improved, with for instance m(q~L) — m{x\) 
being measured with a precision of 4 GeV comparable with the natural widths of the light 
squarks (~ 5 GeV at SPSla) and their mass differences (~ 6 GeV). Further improvement 
in mass measurement precision therefore probably requires that such effects be taken into 
account, for instance in the definition of the x 2 function in Eqn. 3.5. 



4. Conclusions 

A new technique for improving the precision of LHC mass measurements has been outlined 
in which experiment-wise information, for instance invariant mass end-point constraints, 
are combined with event-wise kinematic information such as E™ lss constraints and mea- 
sured four-momenta of visible decay products. For the SPSla model considered here the 
mass measurement precision was shown to improve by ~ 30% for 100 fb _1 of data. SUSY 
models with larger branching ratios for symmetric events containing the necessary decay 
chain, or with more indistinct or poorly-measured kinematic end-points, might be expected 
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to benefit still further from application of the technique. The technique could be extended 
to non-symmetric events where some/all of the SUSY states appear in both legs of the 
event - this strategy for measuring further SUSY masses will be examined in more detail 
in a future paper. Full reconstruction of SUSY particle kinematics with this technique 
could also potentially be useful for other purposes such as measuring the spin-statistics of 
SUSY states. 
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